Deformation of a Trapped Fermi Gas with Unequal Spin Populations 
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The real-space densities of a polarized strongly-interacting two-component Fermi gas of 6 Li atoms 
reveal two low temperature regimes, both with a fully-paired core. At the lowest temperatures, the 
unpolarized core deforms with increasing polarization. Sharp boundaries between the core and the 
excess unpaired atoms are consistent with a phase separation driven by a first-order phase transition. 
In contrast, at higher temperatures the core does not deform but remains unpolarized up to a critical 
polarization. The boundaries are not sharp in this case, indicating a partially-polarized shell between 
the core and the unpaired atoms. The temperature dependence is consistent with a tricritical point 
in the phase diagram. 

PACS numbers: 03.75.Ss, 05.70.Fh, 74.25.Dw 
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The formation of pairs consisting of one spin-up and 
one spin-down electron underlies the phenomenon of su- 
perconductivity. While the populations of the two spin 
components are generally equal in superconductors, an 
imbalance is readily produced in experiments with gases 
of trapped, ultracold fermionic atoms, as was recently 
demonstrated 0,01 • Exotic new states of matter are pre- 
dicted for the unbalanced system that, if realized, may 
have important implications for our understanding of nu- 
clei, compact stars, and quantum chromodynamics. Cal- 
culations show that phase separation between pairs and 
the excess unpaired atoms is one possible outcome for 
a strongly interacting two-component gas 0, 0> IE • 
We previously reported evidence for a phase separation 
in a trapped atomic Fermi gas to a state containing a 
paired central core, with the excess unpaired atoms re- 
siding outside this core 0]. Such a phase separation can 
be detected in the real-space distributions using in- situ 
imaging, where a uniformly paired region produces a min- 
imum in the difference distribution obtained by subtract- 
ing the majority and minority spin densities |2j. 

In our previous work, we noted that the excess un- 
paired atoms reside primarily at the axial poles of the 
highly-elongated trap, while relatively few occupy the 
equatorial shell. As a result, the central minimum in the 
difference images was accompanied by a corresponding 
central dip in the axial density profile obtained by in- 
tegrating the two-dimensional column density along the 
radial coordinate Q|. Several authors have calculated 
spatial distributions for phase separation, assuming both 
a harmonic trappingpotential and the local density ap- 
proximation (lda) jE H, E3, E3, 

It was pointed out 

that, under these assumptions, a uniformly paired core 
would produce a constant axial density difference, rather 
than a central dip 0, Shin et al. have recently 

adopted in-situ imaging, and present their resulting im- 



ages as evidence for phase separation, thou gh in their 
experiment, no such deformation is observed jl5j. 

In this paper, we characterize the properties of the 
phase-separated state using detailed quantitative mea- 
surements of the deformation and reconstructions of the 
three-dimensional (3D) density distributions. Further- 
more, we explore the role that finite temperature plays 
in determining the overall behavior of this system and 
identify two distinct low-temperature regimes. 

Our methods for producing a strongly-interacting, two- 
component Fermi gas of 6 Li atoms have been described 
previously [2J, LLl| - The relative population of two hy- 
perfine states, designated as |1) and |2), is controlled by 
driving radio-frequency transitions between them. Spin 
relaxation is negligible over the duration of the experi- 
ment. A nearly uniform magnetic field is tuned to the 
location of a broad Fcshbach resonance at 834 G G1G1, 
where the two-body scattering length diverges (±oo) pro- 
ducing unitarity-limited strong interactions. The com- 
bined optical and magnetic trapping potential is given by 
U(r,z) = \raoj\7? + U a (l - (^)e~^l^ M) , where 
lub = (27r) 6.5 Hz, w = 26 /mi, w(z) — w [l — (z/z ) 2 }, 
and z — 1.7 mm. Both radial and axial potentials are 
approximately harmonic for sufficiently small r and z. 
The atoms are evaporatively cooled by reducing the trap 
laser intensity until U achieves its final value of 540 nK. 
At this trap depth, the radial and axial trap frequencies 
are ui r = (2ir) 325 Hz and lo z = (27r) 7.2 Hz, respectively. 
The two states, |1) and |2), are sequentially imaged in the 
trap by absorption. The first optical probe pulse breaks 
pairs, causing a small heating that radially broadens the 
second images. We have reduced the delay between probe 
pulses to 27 fj,s, significantly reducing the probe-induced 
broadening previously observed for delays of 215 fis Q. 
Analysis of the images provides measurement of the num- 
ber of atoms in each state, A?i and N2, from which the 
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FIG. 1: In-situ absorption images of a polarized Fermi gas. 
The top figure of each sequence corresponds to the column 
density of state |1), the middle to state |2), and the third to 
the difference of the two. The polarizations are (a) P = 0, 
(b) P = 0.18, (c) P = 0.37, (d) P = 0.60, (e) P = 0.79, and 
(f) P = 0.95. In each sequence, state |2) was imaged first, 
followed by state |1). Slight probe-induced heating can be 
discerned in the images of state |1), where the distribution 
bulges slightly in the radial direction in the region of overlap 
between the two states. The field of view for these images is 
1654 /im by 81 /im. The displayed aspect ratio was reduced 
by a factor of 4.4 for clarity. 
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FIG. 2: Aspect ratio vs. P. The ratio of the axial to the radial 
dimensions, R z /R r , is shown for state |1) by the black circles 
and for state |2) by the red crosses. The radii R r for both 
states are determined by fitting the column density profiles 
to zero-temperature, fermionic Thomas-Fermi distributions. 
The axial distributions are distinctly non-Thomas-Fermi-like, 
so R z is found by a simple linear extrapolation of the column 
density to zero. An aspect ratio of 36 is the expected value 
for a non-interacting gas with anharmonic corrections, in rea- 
sonable agreement with the observations. The uncertainty 
in P is 0.04, which is the standard deviation of polarization 
measurements deliberately prepared as P = 0. There are 
shot-to-shot variations in JVi and a small systematic varia- 
tion towards larger iVi at smaller P. For P < 0.40, iVi = 170 
k ± 40 k, and for P > 0.40, Ni = 135 k ± 25 k, where the 
uncertainty is the standard deviation of the measurements. 
The corresponding average Fermi temperature is Tf ~ 430 
nK, where we define T F = ft(w^ z ) 1 /3( 6i v 1 )i/3/ fcs ^ 



polarization P = (N 1 -N 2 )/{N 1 +N 2 ) is determined. Fit- 
ting the profiles of gases deliberately prepared as P = 
to fermionic nonzero-temperature Thomas-Fermi distri- 
butions, gives fitted temperatures of T < 0.05 Tp, where 
Tp is the Fermi temperature. The actual temperatures 
are expected to be closely related to T . 

Figure 1 shows a series of images corresponding to a 
range of P from to 0.95. The minority spin (|2)) dis- 
tribution becomes markedly less elongated with increas- 
ing P, while its radial size remains approximately the 
same as that for the majority spin (|1)). This deforma- 
tion causes the bunching of unpaired atoms at the axial 
poles, and a lack of them in the equatorial shell. Re- 
markably, the deformation grows continuously with P. 
The central holes in the difference distributions, which 
are approximately equal to the background level for all 
but the highest values of P, indicate that the central core 
is nearly uniformly paired. Figure 2 shows the aspect ra- 
tio for both states. While the majority state aspect ratio 



changes little, that of the minority (representing the core) 
decreases by a factor of 10 when going from completely 
unpolarized (P = 0) to completely polarized (P = 1). 

Figure 3(a) presents a center-line cut of the column 
densities of the majority and minority states, as well as 
their difference, taken along the axial direction. The ax- 
ial density difference profile of this data also exhibits a 
pronounced central dip, as in Ref. 2]. A sharp phase 
boundary between the core and excess fermions is read- 
ily apparent in Fig. 3(a), indicating that a partially po- 
larized shell, if it exists, is extremely thin. Cylindrical 
symmetry of the trap enables reconstruction of the true 
3D density distribution n(r, z) from the column densities 
by use of the inverse Abel transform [2(| • An axial cut of 
the reconstructed 3D density, n(r = 0, z), is shown in Fig. 
3(b). The ratio of the central densities, ni(0, 0)/ri2(0, 0), 
obtained from the reconstructed 3D distributions, are 
plotted vs. P in Fig. 4(a). The central core remains 
unpolarized until at least P ~ 0.9, in contradiction to 
the results reported in Ref. |l5j . where uniform pairing 
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FIG. 3: Column-density profile and 3D density reconstruc- 
tion. The black lines correspond to state |1), the red to state 
|2), and the green to their difference, for P — 0.35 and A?i 
= 175 k. The circles are the result of our theory (see text), 
(a) Center-line (r = 0) axial cut of the column-densities, (b) 
Center-line axial cut of the reconstructed 3D densities. The 
signal to noise of (b) was improved by reflecting and averag- 
ing the column density images about both the r = and the 
z — planes before reconstruction. The densities in (b) are 
slightly diminished due to probe induced radial broadening. 
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FIG. 4: Ratio of the central densities vs. polarization, (a) 
T < 0.05 Tf, corresponding to the data shown in Fig. 2; 
(b) T « 0.2 T F , with average Ni = 500 k. The dotted lines 
correspond to [(1+P)/(1— P)] 1//2 , the expected central density 
ratio for a harmonically confined, non-interacting gas at T = 
0. The increase in ni(0, 0)/ri2(0, 0) in (a) for P > 0.9 may 
be explained by higher temperatures for these data that arise 
from inefficiencies in evaporative cooling at very high P. 



was observed to break down at P « 0.77. 

The observed LDA-violating deformations are surpris- 
ing because the radial size of the distributions is about 10 
times larger than the inverse Fermi momentum Uf, the 
expected magnitude of the correlation length. Nonethe- 
less, several LDA-violating mechanisms, enhanced by 
confinement in a high aspect ratio trapping potential, 
may explain these observations. Gradient terms in the 
Gross-Pitaevskii equation can lead to LDA-violating de- 
formations on the BEC side of resonance, but the magni- 
tude of the calculated effect is much smaller than we ob- 
serve at unitarity |2l|. De Silva and Mueller have shown 
that surface tension between the normal and superfluid 
phases can result in deformations of the minority compo- 
nent that are quite similar to those observed here ■ To 
quantify this, we have developed a theory for the density 
profiles in the gas that allows for a deformed local chem- 
ical potential that is consistent with the macroscopic de- 
formation of the minority density profiles. The results 
of this calculation, the details of which will be published 
elsewhere, is shown in Fig. 3(a) and compares favorably 
with experiment. Within the limitations of our calcula- 
tion, we find no evidence for a deformed Fermi surface 
23] in the experimental data. 

We had previously found that phase separation oc- 
curred only for P > P c , where P c ~ 0.1 2j. For 
P < P c , the observations were consistent with a non- 
phase-separated polarized superfluid. The present data, 
however, exhibits phase separation for arbitrarily small 



P. Since the previous work, we have improved the ef- 
ficiency of the evaporation trajectory, and now obtain 
fitted temperatures that are about half of those previ- 
ously attained. This temperature-dependent behavior 
is consistent with a phase boundary between a phase- 
separated regime and a polarized superfluid (Sarma or 
breached-pair phase) 0, EjJ Ei| at nonzero tempera- 
ture. Such a phase boundary has recently been discussed 
in the context of a tricritical point in the phase diagram 
CHESS EH- To test this hypothesis, we deliberately 
produced higher temperatures by stopping the evapora- 
tion trajectory at a higher trap depth (1.2 fiK), resulting 
in T s» 0.2 Tf- Figure 4(b) shows that the central densi- 
ties in this case remain equal until a critical polarization 
of 0.6 — 0.7 is reached. 

Figure 5 shows absorption images prepared at both 
the lower and higher temperatures. It is readily ap- 
parent from the images that the density distributions 
of the two components of the higher temperature gas 
show no deformations, in contrast to those of the colder 
case. We find that at higher temperatures, the aspect 
ratios of the minority and majority components remain 
equal and constant for all P. This lack of deformation 
is also evident in the axial density distributions, where 
in the case of the colder data, the axial difference dis- 
tribution (Fig. 5(b)) shows the characteristic double- 
peaked structure observed previously Q, while that of 
the warmer cloud (Fig. 5(d )) e xhibits the flat-topped dis- 
tribution reported in Ref. |15| and expected for a paired 
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FIG. 5: In-situ absorption images and integrated profiles. To 
the left are absorption images, while the plots to the right are 
the corresponding axial density distributions. (a),(b): P = 
0.50, Ni = 146 k, with T < 0.05 T F ; (c),(d): P = 0.45, 
Ni = 374 k, with f » 0.2 T F . The central bulge evident in 
the state |1) density in (b) is indicative of strong attraction 
between atoms in the core region. 

core with no deformation 0, . The phase boundary 
is also much sharper for the low temperature data. In 
summary, the higher temperature data support the sug- 
gestion of a temperature dependent transition between 
a low-temperature phase separated state and a higher 
temperature polarized superfluid [2^, |2^| . 

We have reported that pairing with unequal spin pop- 
ulations leads to real-space deformations in a highly elon- 
gated, but still three-dimensional geometry. The sharp 
phase boundaries between the superfluid core and the po- 
larized normal phase are consistent with the usual con- 
vention that phase separation is associated with first- 
order phase transitions. At elevated temperatures, but 
still below the transition to the normal state, deforma- 
tions are absent, and a partially polarized shell is ob- 
served between a uniformly paired core and the fully po- 
larized outer shell. These observations support the exis- 
tence of a tri-critical point in the phase diagram where 
the superfluid-normal transition changes from second- 
order to first-order as the temperature is lowered. 
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